Spin-based all-optical quantum computation with quantum dots: 
understanding and suppressing decoherence 
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We present an all-optical implementation of quantum computation using semiconductor quan- 
tum dots. Quantum memory is represented by the spin of an excess electron stored in each dot. 
Two-qubit gates are realized by switching on trion-trion interactions between different dots. State 
selectivity is achieved via conditional laser excitation exploiting Pauli exclusion principle. Read-out 
is performed via a quantum-jump technique. We analyze the effect on our scheme's performance 
of the main imperfections present in real quantum dots: exciton decay, hole mixing and phonon 
decoherence. We introduce an adiabatic gate procedure that allows one to circumvent these effects, 
and evaluate quantitatively its fidelity. 



I. INTRODUCTION 

The promise of quantum computation is to enable al- 
gorithms which render feasible problems requiring exor- 
bitant resources for their solution on a classical com- 
puter. This has stimulated a large number of propos- 
als for the physical implementation of the elementary 
logical operations building up a general-purpose quan- 
tum computer [J. A central issue is the trade-off be- 
tween efficient coupling to the system, in order to con- 
trol the quantized degrees of freedom, and good isolation 
from the environment, in order to preserve the coher- 
ence of the quantum evolution. Strategies have been 
developed to fight decoherence taking place during the 
computation, ranging from "active" (error-correcting) to 
"passive" (error-avoiding) schemes. Thereby, unwanted 
physical processes (i.e., computational errors) of a general 
kind can be compensated for, either by detecting and cor- 
recting their effect via redundant qubit encoding or 
by decoupling the qubits from the environment dynam- 
ics through algebraic techniques exploiting symmetries in 
the evolution 0]. A less general, more implementation- 
dependent approach is to study the specific decoherence 
channels of a certain physical scheme, and to design gat- 
ing processes that are stable against the relevant types 
of errors (for an example with ion traps, see Ref. 3|). 

In this paper, we adopt the latter point of view, and ap- 
ply it to a recent proposal for all-optical quantum infor- 
mation processing based on charged semiconductor quan- 
tum dots 0]. In this scheme, quantum information is 
stored in the spin of an excess electron in a quantum 
dot (QD), and gating between 2 QDs is performed via 
optical excitation of electron-hole pairs (excitons), which 
in an external electric field acquire a dipole moment al- 
lowing them to interact with each other. In this way, 
the quantum memory coherence time is in the fxs range, 
typical for spin degrees of freedom in semiconductor het- 
erostructures (| , while the two-qubit gating time is in the 
ps range, as dictated by the electrostatic dipole-dipole in- 
teraction. 



Decoherence is important mainly during gate opera- 
tion when excitonic states are created that interact with 
the phonon bath. Recent calculations for the strong field 
limit indicate that the leading dephasing mechanism 
is the coupling to acoustic phonons. In this paper, we 
focus on this coupling mechanism, and we describe a pro- 
cedure that allows one to circumvent, to a large extent, 
the limitations imposed by this decoherence channel on 
the fidelity F of the gate operation. We will evaluate 
the dependence of F on various parameters, including 
temperature, and take into account other sources of im- 
perfection like heavy-/light-hole mixing. 

We choose to neglect all kinds of nonidealities arising 
from limitations, e.g., in QD fabrication and manipula- 
tion techniques. We are well aware that these might yield 
the most significant problems for the implementation of 
our proposal in the immediate future. However, we prefer 
to focus on fundamental quantum-mechanical limitations 
of our physical system rather than on technical problems. 
Once the technological advances have overcome the lat- 
ter, the relevant part will be to find ways to circumvent 
the former. The main purpose of this paper is to develop 
strategies aimed at this. 

The paper is organized as follows: in Sec. II we de- 
scribe the general idea of a two-qubit quantum gate based 
on selective switching of controlled interactions. In Sec. 
Ill we recall the dynamics of charge carriers in a quan- 
tum dot, including external static and oscillating elec- 
tromagnetic fields. In Sec. IV we derive few-level model 
corresponding to the above general scheme and discuss 
some of its limitations. In Sec. V we discuss our two- 
qubit gate and develop its adiabatic version, suitable for 
operation even in realistic scenarios with hole mixing. In 
Sec. VI we propose a hole-mixing tolerant scheme for 
single-qubit operations. In Sec. VII we analyze the ef- 
fect of the interaction with phonons on the performance 
of our adiabatic gates, showing that the gates indeed are 
quite robust also against this kind of imperfection. In 
Sec. VIII we describe how the quantum jump technique 
can be employed for measuring the spin state of a con- 
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fined electron, emphasizing that this can be done even 
for the case of non-zero hole mixing. Our conclusions are 
summarized in Sec. IX. 



II. QUANTUM GATE MODEL 

To execute an arbitrary quantum computation, i.e. to 
control the coherent evolution of a system composed of 
an arbitrary number of qubits, one does not need to re- 
alize physically arbitrary multi-qubit operations. On the 
contrary, just two kinds of elementary operations are suf- 
ficient, out of which all others can be constructed. These 
two elementary gates are the set of rotations of a single 
qubit, and a specific entangling operation on two qubits. 
Among the possible choices for the latter, one which is 
well suited for an implementation with atomic-like sys- 
tems like quantum dots is the phase gate - a transforma- 
tion which rotates by a certain phase just one component 
of logical states: 



|0)|1) 
|1>|0> 
11)11) 



|0)|0) 
|0)|1) 
|1)|0) 

1)11). 



(1) 



When we have -d — tt, this is equivalent, up to single- 
qubit rotations, to a controUed-NOT gate. Ideally, this 
would be accomplished by means of a state-dependent 
interaction of the form 



i7idcal = Ai?,,(t)|l)„(l|®|l),(l| 



(2) 



This describes a situation in which the two-qubit system 
undergoes an energy shift AEat if and only if both qubits 
are in state |1). Imposing the additional condition 



to+T 



AEab{t')dt' = ^ 



(3) 



on the time dependence of the energy shift, Eq. |^ is 
recovered. 



A. Phase gate model: auxiliary interacting states 

An interaction of the form Eq. Q is not straightfor- 
wardly found in nature. Implementing it entails of course 
a certain degree of engineering "natural" interactions, 
i.e., those directly available in a specific physical system. 
This, together with other requirements on the stability 
of the available quantum memory, affects the choice of 
the particular qubit implementation. When it comes to 
systems of confined electrons in solid state systems, such 
as quantum dots, two different choices are natural for 
the logical degree of freedom: either charge excitation 
0, or spin polarization 0. The former provides for a 
strong interaction, leading to comparatively shorter gate 
times but to faster decoherence rates as well; conversely. 



the latter suffers less from the coupling to the environ- 
ment, yielding better stability against memory decoher- 
ence, but bears also a weaker coupling between qubits, 
requiring longer times for gate operation. Aiming at a 
high ratio between coherence time and gate operation 
time leads to conflicting requirements. Reasonable trade- 
offs can be achieved in each case - however, sticking to 
the same degree of freedom for both the memory and the 
two-qubit interaction may be not necessarily the only 
option. For instance, the same effect of the interaction 
Eq. ^ can be obtained by introducing an auxiliary state 
I a;). Let us consider two qubits, labeled by a and b and 
with logical states \a)a,b (a G {0,1}). Each qubit is se- 
lectively coupled to a further state |a:) - namely, only |1) 
can be excited to |a;). This situation is described by the 
following Hamiltonian: 



I/— a, 6 

^Eabit)\x)Jx\®\x),{x\ (4) 



a=0,l,x 
ly—a^b 



As in Ref. the logical states |0) and |1) (the quantum 
memory) can be encoded into long-coherence spin states, 
while the auxiHary states |a;), needed for the gate to be 
performed, can be chosen to be electrostatically interact- 
ing states. State selectivity, required for conditional log- 
ical operations, is accomplished via the state-dependent 
coupling Q(t). The simplest strategy for performing a 
quantum gate exploiting the coupling scheme Eq. (^[1 
would be, e.g., to selectively excite the interacting state 
\x) via a Rabi flop, wait for the desired gate phase to be 
accumulated, and then de-excite. The interaction energy 
shift would then be effective only if both QDs started off 
state |1), as described in Fig. ^ This procedure works 
in the ideal case when the coupling is perfectly state- 
selective as in Eq. (0J). Compared to similar schemes for 
neutral atoms (see e.g. jlfl), it has also the advantage 
that quantum dots, unlike trapped atoms, are not sub- 
ject to back-action on motional dynamics. However, in a 
real situation state selectivity may not be perfectly satis- 
fied, in which case the simple procedure described above 
would not work. We will handle this imperfection below, 
and develop a strategy for overcoming it. But first we 
need a model of quantum dot dynamics that can account 
for Eq. gl). This is the subject of Sect.lml 



B. Gate fidelity 

To evaluate the performance of a quantum gate, one 
needs to compare its desired operation, Eq. in our 
case, with the actual performance of the physical system 
which implements it. The fidelity F represents a quan- 
titative basis for this. To define it, let us start from the 
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FIG. 1: Gate operation via an auxiliary state: in the ideal 
scenario when only one of the logical states can be coupled 
to the interacting state, the interaction leading to the logical 
phase is "switched on" only when both qubits are in the same 
state (here, |1)). 



logical input state 



Ix) = coo |00) + coi |01) + cio |10) + cii 111) = ^ c„ \n) , 

(5) 

which is an arbitrary superposition of all two-qubit com- 
putational basis states. The goal of gate operation is to 
obtain the ideal output 



(6) 



n=0 



This is equivalent to the desired two-qubit transforma- 
tion Eq. JQ): one can be recovered from the other by 
redefining the logical states via single-qubit operations. 
The gate phase "i? turns out to be related to the logical 
phases (j)ap as follows Q: 



1} = <l>00 - (1)01 - 010 + ' 



(7) 



Thus the condition = n simply translates into a condi- 
tion on the (flap's. 

The actual physical situation may involve other exter- 
nal (i.e., non- logical) degrees of freedom, which are not 
perfectly under control. In this case the initial state a 
will rather be a mixture: 



(^ = Ix) (xl «)Poxt(to), 



(8) 



where poxt denotes the density matrix for external de- 
grees of freedom. The operation U realized in the lab 
will in general involve both internal and external degrees 
of freedom in a non-trivial way: therefore the actual out- 
put 



cr' = UaU^ 



(9) 



will no longer be written in a simple factorized form like 
Eq. (jSJ . In order to compare this state with the ideal one 
which would be obtained in the case of perfect operation, 
we define the fidelity 



F = mintroxt (xk' Ix) 

X 



(10) 



The intuitive meaning of this definition is that of a 
worst-case estimate (hence the minimum over the pos- 
sible inputs Ix)) of the gate performance, averaged over 
the available non-logical states not being under control 
(hence the trace over the external degrees of freedom). 
Another option would be to define the fidelity as an av- 
erage over the logical inputs |x) ^ however, we will stick 
to the minimum fidelity which gives a lower bound to the 
average fidelity. 

Taking the most general possible form for the external 
state. 



Pe.t=J2PE\E) {E\ 



(11) 



E 



and assuming that the evolution does not mix different 
logical states, 



W|X)I^) «I]c„|n)®K|S) 



(12) 



n=0 

the fidelity takes the form 
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F = mm V |c„|'|c„|'e^(^"-^'")tre,t(KnT/t), (13) 



m.n— 



which will be relevant for our calculations below. 



III. QUANTUM DOT DYNAMICS 

Quantum dots, due to their discrete density of states, 
are a very promising candidate for the implementation 
of quantum information processing ITU . The brilliant 
idea first proposed by DiVincenzo and Loss to em- 
ploy the spin of an electron confined in a QD as the 
qubit degree of freedom has been developed by the au- 
thors over the years |l2| and is now pursued by many 
groups 0, Q, 01 . Combining QD technology with 
ultra-fast laser pulses now seems to be one of the most 
promising channels for such an implementation scheme 
H Recently the necessary coherence required for 
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such a task, i.e^, Rabi oscillations, has been experimen- 
tally observed There have been also impressive ex- 
perimental achievements in exciting and probing excitons 
in QDs [13. 

The complex many-body dynamics of charge carriers 
in a semiconductor can be considerably simplified when 
considering semiconductor heterostructures like quantum 
wells and dots. The purpose of the present Section is to 
write down explicitly the carrier Hamiltonian for a quan- 
tum dot under these approximations. In the next Section 
this description will be linked to the particular model de- 
scribed by Eq. ^ . Two main approximations 0, |23| 
are understood throughout the following. The first is 
the effective mass approximation, which arises from ap- 
proximating the band dispersion relation around a band 
extremum up to second order in the carrier wavevec- 
tor k. This is valid for small values of k, and allows 
for simplifying Hamiltonians in terms of effective elec- 
tron and hole masses that take into account the under- 
lying many-body dynamics. The other is the envelope 
function approximation |2ll |. which is based on the fol- 
lowing assumptions: (i) the different materials constitut- 
ing the heterostructure are perfectly lattice-matched; (ii) 
the periodic parts of the Bloch functions u/_o(r) are the 
same in the different layers; (iii) the confining potential 
is smoothly varying on the scale of the lattice structure, 
apart possibly from abrupt interfaces. The wavefunction 
can then be expanded as a sum of products of the rapidly 
varying functions u;^o(r) by slowly varying envelope func- 
tions which obey an effective Schrodinger equation in- 
volving the effective masses. We consider QDs in the 
"strong-confinement" regime, in which the typical length 
scale in the growth direction L is of the order of 10 nm 
to 20 nm. Considering QDs in the strong-confinement 
regime means that all relevant energy scales, e.g. charge 
carrier interactions in the QD or electron phonon inter- 
actions, will be small compared to the level spacing of 
the QD, typically of the order of 25 meV for electrons. 



A. Single-particle states under external fields 

Under the above approximations, the carrier Hamilto- 
nian for a quantum dot can be written as I20, 22] 



II - -r -r ^ (14) 

The electron in-plane Hamiltonian iJ| describes the con- 
finement in the direction perpendicular to the QD sym- 
metry axis z, which can be modeled with a parabolic 
potential: 



2nie 



eF • r 



(15) 



where the in-plane coordinate vector is r = {x,y), and 
the electrical field F is taken to be parallel to the xy 
plane. Defining 

re = r+^^ = {re,ee), (16) 



the eigenstates |n,q)^ of iJ| in coordinate representation 
are 



\n,qy 



(17) 



where n = 0,1,... is the principal, q = —n,~n + 
2, . . . , n — 2, n the azimuthal, and rir = {n—\q\)/2 the ra- 
dial quantum number; Cnl (z) are Laguerre polynomials; 



^0 = 



and the eigenenergies are 



h{n + l)u!o. 



(18) 



(19) 



In the growth direction the perpendicular Hamiltonian 
HJ^ is in general a very narrow potential given by the 
quantum well structure and is therefore typically approx- 
imated by a step like potential. In the strong-confinement 
regime, a good approximation is to assume that the sys- 
tem remains in its ground state. Thus the problem effec- 
tively reduces to the in-plane dynamics. The hole Hamil- 
tonian II is of course the same as H'^ || but with hole 
parameters uih and uJh, and opposite charge. 

Due to the strong confinement and spatially symmet- 
ric shapes of the confining potentials in QDs, electronic 
angular momentum states can be defined. They exhibit 
many atomic-like symmetries which have been experi- 
mentally identified 2j\. In contrast to atoms, when con- 
sidering the quantum numbers defining the angular mo- 
mentum of an electron or a hole confined in a QDs, one 
has to take into account the underlying band structure 

0. _ 

Taking spin into account leads to splitting into hole 
subbands. The valence band, built from atomic p-type 
orbitals, contains states carrying an internal (band) an- 
gular momentum m equal to unity. Thus the total an- 
gular momentum is 



j = a + m + l 



(20) 



where a is the spin and 1 the orbital angular momentum. 
Good quantum numbers are the modulus of j and its 
component along the QD symmetry axis z. The single- 
particle states of the valence band with 1 = are classified 
according to the value of (jcr -I- m|, ctj + mz), as follows: 

(3/2, ±3/2): heavy-hole subband; 

(3/2, ±1/2): light-hole subband; 

(1/2, ±1/2): spin-orbit split-off subband. 

For the dynamics considered in this paper, only heavy 
and light holes will matter, the split-off subband being 
energetically far apart. So let us define electron and hole 
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operators for the QD labeled hy v {ly G {a,b}), with 
composite index i — [n, q] and spin a: 



(21) 
(22) 



We can now write the noninteracting part of the carrier 
Hamiltonian for the QD v. 



3/2 



^j,a-f^t,j,cr'^'^d,o-' 



;,CT=±i/2 



j;o-' = -3/2 



(23) 



B. Carrier-carrier interaction 



The electrostatic interaction Hamiltonian is written as 



a, a' 



+ kmmhh (24) 

where the matrix elements of the Coulomb potential 

^2 



y(r-r') 



47re r — r' 



(25) 



are calculated on electron and/or hole wavefunctions, ac- 
cording to subscripts. Here, carrier number conserva- 
tion is assumed, since processes violating this (like e.g. 
Auger recombination and impact ionization) are relevant 
for energies and densities higher than we are consider- 
ing |22|- It should be noted that, in contrast to higher 
dimensional quantum structures, in QDs carrier-carrier 
interactions only induce an energy level renormalization 
without causing scattering or dephasing. 



C. Interaction with a laser field 

Let us consider a laser of amplitude E(t) and central 
frequency impinging on our QD. Under the dipole 
and rotating- wave approximations |25| the corresponding 
Hamiltonian is 



h.C. 



(26) 

where ^fj' is the dipole matrix element between the wave 
functions of an electron with spin a in state i and a hole 
having angular momentum with z component a' in state 
J. The resulting selection rule is that the change in the 
number of electron- hole pairs can only be AN = ±1. 



IV. THREE-LEVEL MODEL 

Let us consider the QD identified by the index v, with 
an excess electron in the conduction-band ground state. 
We label its spin states with 



|0>^ = c|;_o,-i/2 |vac) , 
ID. 



c^,o,i/2 |vac; 



(27) 
(28) 



These are eigenstates of the bare Hamiltonian H^, with 
eigenvalues Cq and respectively, and are not 

affected by the carrier-carrier interaction H^'^. On 
the other hand, in the so-called "trion" - i.e., the 
state obtained from Eas. H27l28|) by creating an exciton, 
having "bare" energy ^ 



/2 



0,-1/2 



more 



than one carrier is present in the QD. Therefore, in 
this case the interaction H^'^ changes the bare state 
''1 +1/2'^!' -1/2^^ ah |vac) into the physical interacting 
state |x, cr^)^, which we will take as our auxiliary state 
for gate operation. Such states were observed and stud- 
ied experimentally in single self-assembled QDs p^ . 

According with the above selection rule, a laser pulse 
can excite at most one exciton. If the laser is tuned on the 
lowest interband excitation energy, then a ground-state 
exciton can be obtained. Due to angular momentum con- 
servation, the hole angular momentum Uh will depend 
upon the laser polarization. For instance, in the case 
of a semiconductor material where heavy holes have the 
lowest energy, and assuming cr'^ circularly polarized light, 
the only hole state that can be excited has — 'i/2. If 
moreover the absolute value of the Rabi frequency, de- 
fined as 



n{t) = 



2/i( 



-1/2, +3/2 



m 



(29) 



is much smaller than the intraband excitation energy, 
|r2| <C uJe,h, then we can neglect the probability of pro- 
moting the electron from the valence band to a higher- 
excited conduction band state. Under these assumptions, 
the interaction Hamiltonian Eq. (|26|l simplifies to 



Vo, 



-l/2"'i/j-,+3/2 



(30) 



If the temperature is sufficiently low with respect to the 
electronic intraband excitation energy, fcsT <C huje, then 
we can neglect also the excited states of the excess elec- 
tron. 

In the subspace defined by the states 
{10},. Jl),. > |a;,-H3/2)^}, the effect of the carrier- 
carrier interaction term H^"^ will be to change the energy 
of the trion states. In particular, an external static 
electric field F applied in the {x, y) plane will mutually 
displace the wavefunctions of the electron and of the hole 
which constitute an exciton, since they have opposite 
charge. In this way, the trion states acquire an electric 
dipole moment. The electrostatic interaction will then 
shift the energy of a trion state where a trion is also 
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present in a neighboring dot. This energy difference, the 
so-called trion-trion shift AEat, will be very important 
for obtaining the state-dependent phase needed for the 
logical gate to be performed. The key ingredient for this 
is a state-selective coupling of the logical states |0) and 
|1) to the auxiliary interacting state \x^') = |x, +3/2). 
In the simplest, ideal case, this state-selectivity can be 
obtained if only one of the logical states is coupled to 
|a;^). The possibility of realizing physically such an 
effect in QDs is offered by the mechanism described in 
the following section. 



static magnetic field. In the adiabatic scheme we are 
going to develop, this will have the effect of suppressing 
unwanted transitions between |0) and |1). This few state 
model will be valid for fl, A, AEab/^, and the Fourier 
width of the pulse t^^ to be much smaller than the QD 
level spacing, so that transitions to excited states are 
negligible. For A{t) = 0, Hcb has the same form as Eq. 

and is thus suitable for quantum gate operation. A 
non-zero A{t) will turn out to be relevant for correcting 
the so-called hole-mixing problem, which is outlined in 
the following section. 



A. Exciton Pauli blocking 

In fact, the Pauli exclusion principle forbids double oc- 
cupancy of any of the electronic states. In particular, if 
the excess electron occupies the state |0), no further elec- 
tron can be promoted from the valence band into that 
state, and thus creation of an exciton by a cr^-polarized 
laser pulse is inhibited (left part of Fig. [2J|. This effect, 
referred to as Pauli blocking, has been experimentally 
verified [23|. On the other hand, if the excess electron 
was in |1), nothing could prevent a second electron from 
being excited to the state |0), thereby creating the trion 
state |a;^) (right part of Fig. ^ . Taking now into ac- 



1-1/2) |l/2) 



1-1/2) |l/2) 



|3/2) 1-3/2) 



^ — 

|3/2) 1-3/2) 



FIG. 2: Pauli-blocking mechanism: a pulse of (j+-polarized 
light can promote an electron from the valence band to the 
conduction-band — 1/2-spin state of a quantum dot only if the 
latter is not occupied, i.e., if the excess electron in the dot is 
in the opposite spin state (right). Otherwise, no excitation 
takes place (left). 



count all of the above approximations and selection rules, 
we can write the following effective Hamiltonian for the 
dynamics of the relevant degrees of freedom in a frame 
rotating at the laser frequency lol- 



u—a,b 



^ ny] 6\1)J1\-A{t)\x^)^^{x^\ 



(31) 



where for generality we have considered a chirped laser, 
i.e. one having a time-dependent detuning A{t) from 



the |1) 



\x 



) transition. We have added a global (i.e., 



independent of the QD label) splitting S between the two 
logical states, which can be realized, e.g., via an external 



B. Hole mixing 

Indeed, Eq. (|31|l does not account for an important 
feature of a real QD system, namely the interaction be- 
tween the hole subbands, described by the Luttinger 
Hamiltonian 28]. In this more accurate description, the 
actual hole eigenstates are no longer the ones described 
above. In particular, the eigenstate of the (heavy) hole 
involved in the dynamics relevant to our study has to in- 
clude a correction from the light-hole state g +1/2 W''^'^) ■ 
Now, a pulse of cr^-polarized light can promote an elec- 
tron from the valence-band state corresponding to such 
light hole into the state |1). This means that the same 
laser we have included in the Hamiltonian Eq. (|31|l has 
a certain probability amplitude to excite an exciton in 
the QD even if the initial excess electron state was |0), 
that is, the laser-coupling selection rules discussed above 
will be weakly violated in a real QD. This effect can be 
included in our simplified three-level model as an addi- 
tional coupling between the states |0) and |a;^), induced 
by the laser with Rabi frequency Q{t) and weighted by 
the effective parameter e whose typical value is at most 
10%. This leads to the model Hamiltonian 



flit) 



(32) 



{\x^)Jl\+e\x^)jO\)+h. 
+ AEab\x^)^{x^ ® ' 
which will be the basis for our simulations. 



V. TWO-QUBIT GATE IMPLEMENTATION 

In this Section we show how the transformation Eq. 
can be realized in practice using quantum dots. We will 
discuss the ideal scenario of perfect Pauli blocking as de- 
scribed in Eq. I|31|) . and then introduce imperfections. 
We will first lay out a strategy to overcome hole mixing, 
based on adiabatically chirped laser pulses. Then, in the 
next Section, we will show that the same strategy allows 
for suppressing the effect of phonon decoherence. 
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A. Ideal gate: e = 



2. Adiabatic passage via dressed states 



In the absence of hole mixing, the laser excitation of 
the trion is perfectly state-selective. In this ideal case, 
gate operation is particularly straightforward. 



1. Direct Rabi excitation 

The simplest quantum gate scheme exploiting the in- 
teraction Eq. (|31|l is based on the following procedure: 

1. selectively excite a trion via a resonant tt Rabi ro- 
tation; 

2. wait a sufhcient time r Trh/AEab for the gate 
phase TT to be accumulated; 

3. de-excite with a second tt pulse to return to the 
logical subspace. 

This is depicted in Fig. El Since the trion-trion splitting 
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FIG. 3: Two-qubit gate via direct Rabi excitation: pulse se- 
quence for exciting and de-exciting the trion state in each dot 
(above); trion-trion population and accumulated interaction 
phase (below). 



can give an interaction time scale of the order of ps, such 
a gate would be pretty fast - however, unfortunately the 
scheme works only in the idealized model of perfect Pauli 
blocking. Therefore we want to develop an alternative ex- 
citation scheme, to be reliable also in the presence of hole 
mixing. This can be achieved by employing an adiabatic 
technique. 



We want to design a process which allows to "switch 
on" the excitonic state for a certain time, and then to re- 
turn to the initial ground state with the highest possible 
probability, by avoiding at the same time spontaneous 
emission. We can achieve this by using an adiabatically 
chirped laser pulse, i.e. one with a slowly changing de- 
tuning from the excitonic transition. Let us start by con- 
sidering our driven two-level system 



17 
2" 



1) (a;| + h.c. 



(33) 



in the strong-coupling regime. Its two eigenstates (the 
so-called dressed states) 



1+) = sin-|l) -f cos-|a;), 
I-) = cos^ll) -sin^lx). 



where 



tan6'=-f7/A (0 < 6i < tt), 
have the energies 

A 1 



(34) 
(35) 

(36) 
(37) 



which are drawn in Fig. 01 Spontaneous emission can 

E/0 




A/0 



FIG. 4: Dressed states in a driven two-level system. 



occur only if the system has a non vanishing probability 
amplitude of finding itself in its excited state. On the 
other hand, for small values of A/f2 the dressed ground 
state contains a significant component of the interact- 
ing state \x). Hence, the following spontaneous-emission 
avoiding excitation procedure can be devised: 

1. the system is prepared in the electronic ground 
state; 



2. the laser excitation starts at a large negative value 
of A/r2, where the lower dressed state tends to |1); 

3. the laser parameters are slowly changed towards 
smaller values of A/il, achieving the transforma- 
tion 



a\0)+P\l) 



|0)+/3 



cos - 1) — sin - 
2 ' ' 2 



(38) 



4. the system is adiabatically driven back to its initial 
state. 

If such a chirped laser pulse is applied to two neighboring 
dots, they will still acquire a state-dependent trion-trion 
phase, due to the admixture from state |a;) that is reached 
starting from state The main constraint is that we 
want to change the Hamiltonian slowly enough such as 
to remain in the lower dressed state of the driven system 
with a high probability. We will show below that it is 
even possible to perform the operation in such a way 
that no phonon mode is excited during the gate, thus 
greatly reducing its sensitivity to temperature. 



B. Hole- mixing-tolerant gate 

In the presence of hole mixing, Pauli blocking does not 
work perfectly. Therefore a 27r pulse for the transition 
|1) \x^) will leave behind some excitonic population, 
causing decoherence. In this case only an adiabatic gate 
operation procedure can ensure that no excitonic pop- 
ulation survive after gate operation. Let us therefore 
analyze it in more detail in the case when e 7^ 0. 



Hole-mixing-tolerant laser excitation 



The Hamiltonian in this case is 



Hi = Ho + S\l)n\ 



f |0> (x| 



h.c. 



(39) 



The level scheme for Hi is drawn in Fig. |S1 It shows an 
avoided crossing of the order of between states |1) and 
la:) like in the two- level well as a much smaller 

one between |0) and \x). However, the latter is found 
at more positive values of the detuning. Therefore, the 
same procedure as outlined above will work also in this 
case, while the adiabatic excitation Eq. I|38(l will be still 
approximately satisfied. 



2. Two-qubit gate via chirped pulse 

Let us now have a closer look to what happens when 
two neighboring dots undergo the above mentioned pulse 
sequence. The two-dot Hamiltonian, including trion- 
trion interaction, is 




FIG. 5: Dressed states in a driven three-level system. A typ- 
ical hole mixing parameter value of e = 0.1 is assumed. 



|a;l),|l3:) 




xl) , \lx) 



FIG. 6: Level scheme of a system of two coupled quantum 
dots including hole mixing. 



The level scheme is depicted in Fig. |H| Although it is of 
course significantly more complicated than the single-dot 
scheme of Fig. |51 the basic feature remains unchanged: 
under the same procedure as in the previous section, 
with the same laser pulse on both dots, different com- 
putational basis states will acquire different trion-trion 
amplitudes. This state selectivity allows one to obtain a 
nontrivial gate phase even in the presence of hole mix- 
ing, still avoiding spontaneous emission as discussed in 
Sec. IV A 21 We carried on a simulation using the follow- 
ing pulse shapes: 



HI + Hi 



^Eab\x)^{x\®\x)f,{x\ 



(40) 



n{t) = n^e 



A(t) = Aoo 1 



'(t/rfi)" 



(41) 
(42) 
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Results of the simulation with A_Eaf, = 2 meV, UAoo = 
3 meV, tq = 10 ps, ta — 8.72 ps, S = 0.5 meV and 
hilo = 3 meV, are reported in Fig. The gate phase 
§ — n in obtained in a longer time than in the simple 
Rabi-flopping scheme (Sec. IV A since the procedure 
now has to be adiabatic to avoid excitations to decay- 
ing states. Indeed, the population left in the unwanted 
excitonic states remains below 10^^. 



laser intensity 




/ detuning ~ \ 



I .( 



I 



I 



\ 



I Jl \ 
I ^/ \ 
/ / \ 

populations/' \ 



Ti; 



-10 10 

time (ps) 

FIG. 7: Hole-mixing-tolerant two-qubit gate operation via a 
chirped laser pulse: pulse shapes (above); interaction phase 
and trion-trion populations (below). The dashed line depicts 
the population in the \xx) state obtained starting from |11); 
the flat solid line the ones obtained starting from |00), |01) or 
1 10). 



VI. SINGLE-QUBIT OPERATIONS 

Contrary to atomic QC implementation schemes, im- 
plementing the single-qubit gate employing the spin state 
of an electron in a QD is of greater difficulty than imple- 
menting the two-qubit gate p9|. A natural candidate 
for an optical implementation of the spin rotation would 
be employing a two-photon Raman process involving an 
intermediate hole state. However, such a scheme is not 
possible since the lowest lying (long lived) hole states 
have Mj — ±3/2 symmetry, and there is, however, no 
strong dipole allowed two-photon coupling connecting the 
qubit states. This problem can be overcome by applying 
a transverse magnetic field, mixing the spin states and 
addressing the Zeeman split states in a frequency selec- 
tive way, which however severely limits the Rabi frequen- 
cies and thus the gating time An alternative is to 



make a Raman process via the lowest light hole states 
Mj = ±1/2 which however, being excited hole states, 
suffer from significant decoherence. 

Employing II- VI semiconductors grown QDs avoids the 
above mentioned decoherence problems, since in these 
QDs the strain can shift the energy of the light holes to 
become the energetically lowest hole states 30] . Exploit- 
ing such QDs the single-qubit gate can be performed by 
the following pulse sequence: (i) a linearly polarized laser 
pulse couples the light-hole subband and the bottom of 
the conduction band (see Fig.|HJ|. This linearly polarized 



1-1/2) |l/2) 



2 (1 



1-1/2) |l/2) 



|l/2) 1-1/2) 



|l/2) 1-1/2) 



FIG. 8: Single-qubit operation. 



pulse, which can be described as an equal-weighted super- 
position of (7+ and cr_ polarized light, attempts to create 
a trion-trion state in which both ground state electronic 
spin states are occupied but due to Pauli principle, a tt- 
pulse of such a light will promote an electron only into 
the unoccupied qubit state - processes (1) in Fig. |H1 (ii) 
A further 7r-pulse of light now recombines the hole 
state with the original excess electron Fig.|S| In order to 
perform this recombination, i.e., changing the total angu- 
lar momentum j by one unit and keeping its component 
along the QD symmetry axis z fixed, the laser pulse is 
shined in in-plane direction. This is achieved by employ- 
ing a wave guide based scheme ..31j . The Hamiltonian 
describing the above process is given by: 



H 



iq 



?^^(l-/^>(0|+e|xj)(0| 
+ \xi^){l\+e\xi^){l\+h.c.) 

+ a^(|xj)(0| + |xi^)(l|+h.c.; 
~A{t){\xi^){xi^\ + \xi^){xi^\) , 



where \xi^) = \xi,+l/2) ,\xi^) = |a:j,— 1/2) the I index 
notation is used to indicate that the excitonic states are 
defined in terms of light holes, rather than of heavy holes 
as before. The electron-electron interaction appearing in 
the intermediate state has been absorbed by redefining 
the detuning. One should note that the heavy-light hole 
mixing has an effect only for the laser pulse directed in 
the growth direction, z, the in-plane directed laser pulse 
does not induce transitions for the heavy hole part of the 
mixed wave function due to symmetry considerations. By 
properly adjusting the duration and the phase of the laser 
pulses, i.e., adiabatically eliminating the excitonic states, 
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the following effective Hamiltonian is obtained: 

<^ = ?^^^|l>(0|+h.c., (43) 



where fie// = y/(^h^h + £^W)/^ the effective Rabi 
frequency for the coherent rotation between the two logi- 
cal states. Thus utilizing strain inverted heavy-light hole 
II- VI semiconductors grown material, on which the adi- 
abatic two qubit gate procedure can be adopted, single 
qubit optical gating can be induced on a pico-second time 
scale. 



VII. PHONONS AS A SOURCE OF 
DECOHERENCE 

A. Interaction of quantum dots with phonons 

The fidelity of the proposed qubit turns out to be quite 
good, and is determined by several mechanisms. Since 
the expected result is quite small, we can consider differ- 
ent sources of infidelity separately. 

Due to the potentials confining electrons in all spa- 
tial dimensions which lead to a discrete density of states 
QDs are often referred to as "artificial atoms" . The ma- 
jor difference with respect to atoms is the coupling of the 
electrons to underlying lattice degrees of freedom, which 
lead to relatively faster decoherence times, on the order 
of tens of picoseconds [s^l ■ In what follows we are devel- 
oping a simple model of a QD interacting with a thermal 
bath of phonons. 

In a large system of volume V the displacement field is 
linear in creation and annihilation operators of phonons, 
bjq and &jq respectively. Every phonon mode is charac- 
terized by its polarization j, momentum q and frequency 
LUj (q) . The number of phonon degrees of freedom is lim- 
ited by the total number of atoms forming the lattice. 
Hence, there is a maximum frequency, usually defined 
via Debye temperature 8, so that tfj(q) < Q (hereafter 
we use units such that ?i = fcs = 1). Normally O is very 
high (hundreds of K). 

The minimal model describing the interaction of a 
charged quantum dot with the phonon field is given by 
the following Hamiltonian: 



ifph = i^QD - MoE(ct/it + he) + J2u;,{q)b]^b 



'jq Vint 



(44) 

The quantum dot is assumed to be in a time-dependent 
laser field E. Mq is the dipole moment, and the Hamil- 
tonian of an isolated QD is given by 



(45) 



where c and h are the operators annihilating a ground- 
state electron and hole, respectively. Their spin depen- 
dence is not relevant here and will not be considered in 
what follows. As above, the quantities and are the 



energies of the single particle electron and hole states, 
respectively. The phonon coupling Hamiltonian is 



E 



[(5|qC^C - g'l^h^h)bjq + C.C.] , 



(46) 



where g^^ and g^^ are the coupling constants. The 
dynamics only couples the states \g) = |vac) and 
|e) = c^/i^jvac). The matrix elements in this basis are: 
{e\c^c\e) = {e\h^h\e) = 1 and {e\c^h''\g) = {g\ch\e) = 1. 
These states are basically the logical and auxiliary states 
of the three- level model discussed in Sect. IIVI without the 
excess electron, namely, |0) = _i/2\9)' 1^) ~ "-o 1/2 Iff)' 
and \x,a) = clje). 

Rewriting the laser field in the form: E = 
E (t) cos{ujot) , where luq is the laser frequency and E{t) 
is the slow envelope, in the rotating wave approximation 
we obtain: 



ph 



-A + EA,q(6,q + 6]J|e)(e 



n 



+ 2 (|e)(5l + |g)(e|) +^c.,(q)6]q&,q, (47) 
j:q 

where A is the laser frequency detuning, fl = |MoE(t)| is 
the Rabi frequency, and the phonon coupling strength is 
given by Ajq = 9jq~9j'q- This is identical, apart from the 
cavity terms, to the Imamoglu-Wilson-Rae Hamiltonian 
Hi (also used in Ref. 0). 

As we shall see, for slow processes considered below 
the interaction of a quantum dot with optical phonons 
is practically irrelevant. Indeed, optical phonons have a 
gap: u!opti<i) — >■ ojopt 9^ at g ^ 0, which is large and 
hence optical phonons can always be adiabatically elim- 
inated from the low frequency dynamics of the quantum 
dot. If taken into account, optical phonons contribu- 
tion only renormalizes some quantities in the Hamilto- 
nian H47() . This is of course not really important for us, 
since we consider all the constants in the Hamiltonian 
as being phenomcnological (essentially taken from ex- 
perimental data). On the contrary, longitudinal acoustic 
phonons (LA phonons) have no gap: w(q) — uq at q ^ 
(u is the velocity of sound) and can effectively interact 
with low-frequency degrees of freedom of the Hamilto- 
nian (|47|l . Hereafter we will only consider LA phonons 
and omit the index j everywhere. 

The Hamiltonian (|47|l formally coincides with the 
Hamiltonian of the dissipative spin-boson model. The 
latter describes the interaction of a two-level system 
(spin) with the bath of harmonic oscillators (bosons). As 
described in 34] , the important properties of the interac- 
tion of the quantum dot with the phonons are contained 
in the integrated quantity 



(48) 



One of the most important properties of the spectral 
function J is its frequency dependence at small lu: 
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J{lo) ^ uj". The different values of the exponent s dis- 
tinguish between the cases of ohmic (s = 1), sub- (s < 1) 
and superohmic (s > 1) couphngs. 

The phonons are coupled to the charge distribution in 
a quantum dot by means of either deformation or piezo- 
electric coupling potentials. The calculation of the spec- 
tral function J requires a specific microscopic model. In 
the simplest case of a quantum dot characterized by a 
harmonic confinement potential the calculation is pretty 
easy, even if the QD is placed in external electric field 
(see Appendix IXjl . The results of the calculation can be 
summarized as follows: both in the case of deformation 
and piezoelectric coupling the spectral function is super- 
ohmic, with s — 3 and s — 5, respectively. In both of the 
cases the spectral function can be approximately written 
as 



J{ll>) uj" exp(-u;^/tjf ) 



(49) 



with a cut-off at uji ^ u/l, where / is the size of the 
quantum dot. This frequency is nothing else but the 
inverse phonon flight time through the quantum dot. The 
electric field (of reasonable intensity) does not change the 
exponent s of the spectral function. 



B. Adiabatic Hamiltonian: Dressed states 

Instead of considering the "bare" states |e) and \g) it 
is convenient to switch to the "adiabatic basis": let us 
diagonalize first the quantum dot part of the Hamiltonian 
(|47|l . The eigenstates (in terms of the bare states |1) and 
which in the present discussion are replaced by \g) 
and |e), respectively) and energies are given in Sec. IV A 2l 

The full interacting Hamiltonian (|47|l can be rewritten 
in the new basis and split into two parts: ifph — Hd + H', 
where 



H, = 



E+ + J2 ^i<lKb^ + cos' 2^^ibl + 6q)] 1+) (- 
q 

E_ + ^c^(q)5^5q - sin' ^Aq(6^ + b^)] \~){- 



and 



H' 



sin 6* 



(50) 



E-^<i(^l + M(l+)(-l + l->(+l), (51) 



where tan(6') = —il/A so that 6* is a time dependent 
quantity. The roles of the two Hamiltonians Hd and H' 
are very different and will be considered separately below. 

The interaction with the phonons can be understood 
in two ways. One of the remarkable non-perturbative 
simplifications can be used due to the fact that we are 
only dealing with the superohmic coupling case. For a 
slow process (the Hamiltonian parameters change on a 



time scale tcui ^ 1) there is a way to adiabatically elim- 
inate most of the "fast" phonon degrees of freedom with 
a;(q) ^ uji without relying on the pertubation ex- 



< 



pansion in phonon couplings A^. The resulting effective 
Hamiltonian has the same form as Eqs. (|50I51I) . but with 
the summation restricted only to the phonon modes with 
i-u{q) ^ cj* and renormalized values of the Rabi frequency 
Cl and the detuning A (see the Appendix IB)i. In what fol- 
lows we will use both representations on the same footing 
and make no distinction between the bare and the renor- 
malized quantities whenever it is not relevant. 



C. Diagonal and off-diagonal channels; 
Landau- Zener theory 

Our qubit proposal relies on the fact, that the quantum 
dot stays in the same adiabatic "dressed" state 
under slow variations of external parameters {Q or A). 
Therefore, transition between the adiabatic states H37() 
are a source of infidelity. Realistically every gate opera- 
tion is performed with a finite speed and thus undesired 
transitions between the dressed states are always possi- 
ble. In the absence of phonons the transition probability 
is given by the Landau-Zener theory. In its simplest ver- 
sion, i.e. for a case of linear detuning sweep A = At 
around the resonance value A = 0, the measure of infi- 
delity is given by the probability 35] 



P± = exp(-7rf7V4A). 



(52) 



By requesting this quantity to be small we establish our 
adiabatic condition 



(53) 



The condition has a simple physical meaning: the reso- 
nance is observed approximately when we have A ~ 17, 
so T ~ fi/A is nothing else but the characteristic time 
of the detuning sweep. Then, the adiabatic condition 
naturally implies that we have fir ^ 1. 

The effects of the interaction with phonons on this de- 
coherence channel do not change this result much. As 
discussed in the Appendix ^ both in the case of per- 
turbation theory and in the adiabatic approximation the 
phonon interaction only renormalizes the Rabi frequency 
[see Eq. (IB4|l ]. Therefore, the same expression (|52|) . but 
with the renormalized value of D, instead of f2, holds even 
if the phonon coupling is strong. 

Now let us consider the interaction with the phonons 
in some more detail. The phonon-assisted transitions be- 
tween the dressed states |±) are possible and are de- 
scribed by the Hamiltonian (|51|l . Since most of the 
transitions occurs close to the resonance A = 0, when 
the characteristic energy difference between the adiabatic 
levels is ~ ri, and the inter-level transition probability 
can be estimated using the Hamiltonian H' and the Fermi 
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Golden Rule 



|T±(w)pJ(w-17)(iu;, 



(54) 



where 

T± i^) = \ j exp(jwt) sm[e{t)] (55) 

is the Fourier component of the "couphng potential" in 
Eq. (|51|l . The latter can be easily estimated in the adi- 
abatic limit by using the linear approximation for the 
detuning close to the resonance point (as above): 



exp(-t j/tJm) 

{uJUJr, 



U/2 



(56) 



where we have LOm = = A/f2, and lo ^ LOm- The 
frequency LOm is the high-frequency cutoff imposed by the 
speed of the frequency detuning sweep. Substituting the 
above expression into Eq. (|54|l and integrating over cj, we 
find the following estimation for the interstate transition 
probability: 



P± 



n 



(57) 



where a ~ 1 is a numerical factor. This quantity charac- 
terizes the probability of unwanted processes and hence 
is a measure of infidelity. The result is exponentially 
small if we have ^ f^, which is nothing else but our 
adiabaticity condition H53I) . 

At finite temperatures there is an additional mecha- 
nism for decoherence: a quantum dot can interact with 
the phonon field and absorb a thermal phonon. Ac- 
cordingly, the quantum dot acquires energy and is trans- 
formed into the excited dressed state. Nevertheless, the 
process can be easily suppressed by operating in the 
regime of small temperatures T <^ Q. In this case the 
transition probability is characterized by an additional 
small factor exp(— il/T) <C 1 and can be disregarded. 

Both probabilities H52|l and (|57|) have similar structure 
and are exponentially small for adiabatic processes (|53|l . 
In later sections we will find that the "diagonal" terms 
Hd in the quantum dot Hamiltonian, though not chang- 
ing the adiabatic states of the quantum dot, lead never- 
theless to excitations of acoustical phonons and thus to 
a certain infidelity. The results do not contain exponen- 
tially small factors and hence the off-diagonal terms in 
the Hamiltonian considered here can be neglected in our 
further fidelity calculations. 



D. 



General expression for the fidelity in the 
presence of phonons 



The major source of infidelity is the excitation of acous- 
tical phonons without change of the quantum dot state, 
i.e., pure dephasing. To develop a formal approach to the 



fidelity calculation, we consider the evolution of a quan- 
tum dot coupled to a heat bath of phonons. Assume that 
at < = — cx) the system starts from the state which is a 
direct product of pure quantum dot state and a thermal 
state of phonon field at a temperature T. As discussed in 
the previous section, the dynamics of the quantum dot 
can be well described (up to a few exponentially small 
amplitudes) by the diagonal Hamiltonian (|50|l . Using 
this simplification, we can rewrite the density matrix of 
the quantum dot subsystem at i = oo as 



p^p^(UlTTUp)\a){(3\, 



(58) 



where \a) is the quantum dot state, (...) is the average 
over the initial phonon state, T is the time ordering sign, 
and the (diagonal in the dressed state basis) evolution 
operators are given by 

rcxp{-z f dt[/a(t)^Afe(6fe + 6^) + S„]}, 

k 

(59) 

with /+ cos2(6l(t)/2) and /_ = - sm^{e{t)/2). The 
time ordering is not convenient and can be removed by 
transforming the evolution operator into 

r[/„ = exp|-z / dt\u{t)Y,Mbk+bi)+EJ+i^A, 

(60) 

where 

dt dt'f^{t)f^{t')Y,>^l&^^^k{t-t') (61) 

-oo ^oo ^ 

is the phase originating from the non-commutativity of hk 
operators at different times. Combining these results to- 
gether we obtain the following expression for the fidelity 

matrix T„/3 = {UlTTUp): 

Tap = ( exp| -I / dtaapit) ^ \k{bk + b\) + icfj^p > 

(62) 

where a„/3 = fa- fp, and 

/OO 
dt[Ef3{t)~Ea{t)] (63) 
-oo 

+ JJ dtdt'fa{t)fp{t')Y,>l sincukit-t'). 

k 

The fidelity matrix T can be used to rewrite Eq. H58|l 
in a more convenient form. Consider an arbitrary pure 
state of a quantum dot '5 = Ca\a), which evolves into 
the density matrix 



Pat3 = C*^Cf3Taf3{Xk)- 



(64) 



Without the interaction with phonons the evolution is 
characterized by the T-matrix with all Afc set to zero. 
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Therefore, one can characterize the phonon interaction 
by the degree of infidehty, defined as 



f=l-F= 



af3 



(65) 

This is a standard problem of hnear optimization, whose 
solution can be done in the general form: the infidelity is 
given by the largest eigenvalue of the matrix T(A) — r(0). 

In what follows we will estimate the value of the infi- 
delity / for single qubit operations and for the quantum 
gate realization proposed above. 



E. Fidelity of Rabi rotation 

We consider first the simplest case and calculate the 
fidelity of single qubit operations, such as a reversible in- 
ternal rotation. To be specific, we calculate the fidelity 
of an adiabatic sweep of the detuning A around its res- 
onant value A = 0, while keeping the Rabi frequency f2 
constant. Using our fidelity definition from Eq. (|65f) . we 
find that we have / = 1 — cxp(— F), with 



2iVH], 



(66) 



where N{uj) = [exp(ti;/T) — 1] ^ is the phonon occupation 
number, and 



If the coupling with the phonons is weak, i.e. F <C 1, 
then the infidelity coincides with F and is only a power 
law small (compare Eqs. (|69|l and (|70(l with the expo- 
nentially small results 1)52(1 and (|57|l of our off-diagonal 
Hamiltonians discussion). At zero temperature the con- 
ditions F <C 1 and the perturbation theory expansion 
parameter (|B5p are the same. 

Eqs. H69(l and (|70|l are obtained assuming that Um <C 
u!i. In this case the results are not confined to the per- 
turbation theory limit (see the Appendix IbI for more de- 
tails about the adiabatic elimination of high frequency 
phonons) . In the other limiting case the infidelity is given 
by the same Eqs. (|69|1 and (|70l) but after the substitution 



Of course, the condition 



> 



LUi breaks the 



adiabatic separation of the slow and fast phonon degrees 
of freedom. Therefore, the results for the fidelity in this 
regime can only be valid if the infidelity F is small. 

Eq. H68|l is derived in such a way, that its validity is 
not confined solely to the analysis of acoustical phonons. 
In fact it also allows one to understand how the contri- 
bution of the higher frequency degrees of freedom (such 
as optical phonons) can be ruled out. Indeed, optical 
phonons are characterized by the minimum frequency coq 
(the optical gap), so that J{uj) — for all uj < luq. Sub- 
stituting this definition into Eq. H68|) and integrating in 
the adiabatic limit <C loq we find 



F exp( ), 

LUq LUm 



(71) 



dtex.]}{—iu!t) cos{0). 



(67) 



In order to analyze the QD dynamics and compare the 
results with the discussion of the off-diagonal processes 
we use the same sort of linear approximation for the time- 
dependent detuning A = At close to the resonance. A 
simple calculation gives 



where Jq = J(wo)- This is once again an exponentially 
small result (compare with Eqs. 1)69(1 and 1(70(1 with a 
clear physical meaning: a slow process occurring on a 
time scale uj^^ can not excite high frequency lattice vi- 
brations if LUq ^ UJm- 



Rabi oscillations 



1 



A2 ' 



2N{Lo)]duj, (68) 



where we have uJm = A/fl fl (see the discussion 
above), and Ki{x) is the Bessel function of the second 
kind. The main contribution to the integral originates 
from the range of frequencies uj < uim (assuming, of 
course, Um ^ ^/)- Therefore, at small temperatures 
T LUm we can neglect the thermal occupation of the 
phonon states and find that 



F 



J{ujjn) 



UJr, 



(69) 



In the opposite limit, i.e. when T 3> uj^, the phonon 
numbers can be approximated as N{u!) « T/uj and the 
integration yields 



J{uJrn) T 



(70) 



Another revealing example of phonon interaction ef- 
fects is the damping of Rabi oscillations. Consider the 
case of exact resonance (A = 0) and a quantum dot start- 
ing at t = in the state \g). In the dressed state picture 
this corresponds to 



l5) = ^(l+> + |-)). 



(72) 



As time progresses, the state changes and the probability 
to find the quantum dot in the state \g) can be found 
using the density matrix from Eq. 1(64(1 



Pr, 



-[T++ +T__ +23fir+_]. 



(73) 



In our diagonal approximation, at i ^ oo, this is equiv- 
alent to 



Pg ^ - [1 + cos(nt + 4)+- 



(74) 
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with given by the first two terms in Eq. (|63|l . This 
means that the diagonal Hamiltonian H5U|1 describes un- 
damped Rabi oscillations at the frequency fi. Within the 
adiabatic approximation Q <^ tui one can integrate out 
the high frequency phonons (see Appendix ^ and ob- 
serve that in the first approximation the effects of phonon 
interaction show up in the renormalization of the Rabi 
oscillation frequency ^ i7, as given by Eqs. ljB4|) . (|B6|I . 

The gradual damping of the Rabi oscillations origi- 
nates from the off-diagonal Hamiltonian H' (|51l) . In con- 
trast to our previous discussion of the internal qubit ro- 
tation, in this case there is a finite probability to find the 
quantum dot in its excited state |-|-). This means that 
now the processes leading to emission of phonons become 
possible. The transition rate T can be calculated using 
the Fermi Golden Rule 



100 



(75) 



so that 



Pa~2^^ + cos(^^i + <?i+-) exp(-ri)]. (76) 

Since realistically it is <C loi, we have the ratio 
J(ri)/r2 <C 1 and thus the quantum-dot oscillations are 
only weakly damped. 

Altogether this lets us conclude that in the presence of 
phonons a quantum dot in an external laser field under- 
goes weakly damped Rabi oscillations, characterized by 
the renormalized frequency and the damping rate deter- 
mined by the spectral function J(0). 



G. Quantum Gate fidelity 

The ultimate goal of our calculations is the fidelity of 
a quantum gate. The Hamiltonian of a couple of inter- 
acting QDs can be represented as follows 



n 



J2i\e)A9\ + \9)Ae\) + Y.^(^KK 



+AEab\e)a{e\ «) |e)h(e|. 



(77) 



where the index ~ a,b labels the quantum dots, 
Xi, = ±d/2 is the position of the dots, and the last term 
represents the trion-trion shift AEab- 

The trion-trion shift is a crucial element of our quan- 
tum gate proposal. Its presence introduces the condi- 
tional dynamics. The timing of the adiabatic process 
should be designed in such a way that the trion-trion 
shift induces the required phase shift of the state |e)a|e)b. 

The whole analysis of a single quantum dot operation 
can be easily generalized to the quantum gate case. As 
discussed above, we first transform into the dressed state 
basis and then select only "diagonal" terms from the in- 
teraction with the phonons. Then one can calculate the 
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FIG. 9: The fidelity of the quantum gate as a function of 
the inter-quantum dot separation d, for difi'erent values of the 
temperature T: K (solid line), 5 K (long-dashed line), 15 K 
(dashed line). 



fidelity matrix T (now 4x4) and find the fidelity from 
Eq. (|65|l . This program can be completed numerically. 

Let us consider first the two simple limiting cases. In 
the simplest case of a small trion-trion shift AEab <C ^, 
the two quantum dots can be considered separately and 
the trion-trion shift can be accounted for as a perturba- 
tion. Then, following the steps of our single quantum dot 
calculation, we find that the infidelity is / = 1 — exp(— F), 
with 

F = J J{Lj)dLJC0S^{u;d/u)[l + 2N{u;)]\a{uj)f, (78) 

where a(a;) is again given by Eq. 1)67(1 . As expected, in 
this limit the result is practically the same as we had for 
a single qubit operation - see Eq. (|^ . 

The dependence of the gate fidelity (|75|l on the quan- 
tum dots separation d is very weak (the effective value of 
the cos-function under the integral sign is anything be- 
tween 1/2 and 1 for large and small values of d, respec- 
tively). This means that the obtained result is practically 
insensitive to the assumptions regarding the degree of co- 
herence between the phonon modes around each of the 
dots. Indeed, Eq. (|78|) implies that both of the quantum 
dots interact with the same phonon bath. This corre- 
sponds to a case of having the dots interacting with the 
same set of bulk modes. Another possibility is that the 
inter-quantum dot separation is smaller than the phonon 
mean free path. In the case of separate phonon baths 
(which means the quantum dots electronic degrees of 
freedom interact with independent phonon modes, or the 
separation between the dots exceeds the mean free path 
of phonons) , the quantum dots become completely sepa- 
rate and the infidelity is given by its single-qubit expres- 
sion H66|) . Both expressions are indistinguishable within 
a factor 1. The weak dependence of the fidelity on the 
interdot separation at not-too-high T can be seen on the 
Fig. 13 obtained by our exact numerical calculation. 

The calculation in the other limiting case AEab ^ ^ 
is very similar. In this case the quantum gate has two 
avoided crossings instead of one and their contributions 
add independently. Since the width of the resonance is 
in both cases fl, the contributions of independent reso- 
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n (meV)|/(r = 0.1 K)| /(r = 4 K) | /(T = 20 K) 
0.3 I 1 2.2^1.9 X lO'^ l s. 9 H- 4.5 X 10'^ 
0.06 9.7^4.6x10-^4.8^2.3x10"'' 

0.02 7.1 5.1 X 10^^ 5.3 ~ 2.5 X lO-"^ 



TABLE I: The infidelity f(T) as a function of the tempera- 
ture T for different values of the Rabi frequency Q. The dot 
size is 20 nm. The trion-trion shift AEab varies from 0.1 meV 
to 3 meV. 



!^ (meV)|/(r = 0.1 K)| /(T = 4 K) |/(r = 20K) 
0.3 I | l.4^0.92 X 10~^ | 0.15-0.1 
0.06 1.1 X 10"* 1.2 X 10"^ 
0.02 N/A N/A N/A 



TABLE III: The infidelity /(T) as a function of the temper- 
ature T for different values of the Rabi frequency Q. The dot 
size is 10 nm. The trion-trion shift AEab varies from 0.1 meV 
to 3 meV. 



Q. (meV) 


/(T = 0.1 K) 


f{T = 4 K) 


/(T = 20 K) 


0.3 





5.6^4.1 X 10"^ 


0.034 ^ 0.020 


0.06 





5.0^4.6 X 10-'^ 


2.5 ^ 2.3 X 10"^ 


0.02 





1.5 ^ 1.4 X 10"^ 


7.2^6.9 X 10"'^ 



given by 



TABLE IL The infidelity f{T) as a function of the tempera- 
ture T for different values of the Rabi frequency Q. The dot 
size is 15 nm. The trion-trion shift AEab varies from 0.1 meV 
to 3 meV. 



nances add separately and are again of the same order of 
magnitude as Eq. H69(l and H7Q|I . The dependence of the 
fidelity on the interdot separation is very weak again. 

To quantify the discussion above we performed the ex- 
act numerical calculation of the fidelity matrix Ht)2|l using 
a certain shape of the detuning sweep A(t). The fidelity 
as a function of the interdot separation and of the tem- 
perature for a specific value of the trion-trion shift is 
plotted in Fig.lO. The figure nicely shows both temper- 
ature regimes (IS^ and {TD)), as well as the weak depen- 
dence of the result on the separation between the dots. 
We also performed a few calculations for quantum dots 
of different sizes (effectively varying the cutoff parame- 
ter uji). The results of the calculations are summarized 
in the Tables The interparticle separation is 5 

nm, the pulse duration is 1 ps. 



H. Discussion 

In this section we performed a systematic study of var- 
ious decoherence mechanism associated with the interac- 
tion of the quantum dots with phonons. The results of 
our study can be summarized as follows. 

The details of the interaction with phonons can be 
"compressed" into the spectral function J(w). It is char- 
acterized by the strength of the coupling, the frequency 
dependence at small uj and the value of the high frequen- 
cies cutoff uji ~ u/l, where u is the velocity of sound and 
I is the size of the quantum dot. 

The infidelity / turns out to be quite good for all the 
realistic situations we considered. This means that the 
perturbation theory is well applicable and, in the limit 
of small temperatures, the infidelity / = 1 — exp{— F} is 



J{luc) 



(79) 



where lOc = min(ct;;,a;m) and ~ A/ 57 is the inverse 
characteristic time of the gate operation. In the higher 
temperature limit (T ^ ojc) the infidelity scales linearly 
with the temperature 



(80) 



These expressions are the central results of the section. 
They can be applied to estimate the fidelities of both 
single qubit operations and of the quantum gates (in the 
latter case there is also a weak dependence on the sep- 
aration between the quantum dots). The accuracy of 
the simple estimations along the lines of Eqs. H79(l and 
(|80|l was checked by exact numerical calculations of the 
fidelity in a wide parameter range. The contributions of 
the piezoelectric coupling are numerically smaller by 2-3 
orders of magnitude in all our calculations. 



VIII. STATE READ-OUT BY QUANTUM 
JUMPS 

A necessary requirement for a quantum information 
processing implementation scheme is the ability to per- 
form an accurate measurement of a single qubit. Imple- 
mentation of a highly efficient solid-state measurement 
scheme designed to measure the spin or charge of single 
electron is a highly difficult task jSg- 

Monitoring the fluorescence from a single quantum dot 
(QD) has been suggested as a mean to measure single 
scattering events within QDs [s^ , as well as a means for 
final read out of the spin state for the purpose of quan- 
tum computation Recently, it has been verified ex- 
perimentally that the spin state of an electron residing 
in a QD can be read using circular pumped polarized 
light |38j| . In this section we describe how it is possi- 
ble to devise an optical read out scheme based on the 
idea of the Pauli-blocking in QDs even in the presence of 
heavy/light hole mixing. We describe two situations: an 
ideal case with no heavy/light hole mixing and the re- 
alistic case which includes mixing. In the ideal case the 
time of measurement, i.e. the time in which one can still 
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extract the information regarding the spin state of the 
confined electron, is hmited only by the spin decoherence 
time whereas in the case of mixing the measurement time 
is also limited by the typical time for a spin flip induced 
by the excitation process. 

The system we have in mind is described by Fig. 1101 
It is the single-QD counterpart of Eq. (|32|l . taking into 
account also the decay rates from the excited level, which 
we call \x) here for simplicity. We ignore the detunings 
6 and A, which are not relevant here, as well as the de- 
cay rate Kojfrom state |0) to |1), since the typical time 
scale for it [23 is much larger than the typical time scale 
in which the spin directional information is lost due to 
the laser mediated spin flip which is now the bottleneck 
process limiting the measurement time. 




|0) 



ll) 



FIG. 10: The lambda configuration one has to consider to 
include hole mixing. 



A. Ideal case, i.e. no mixing 

Let us first consider the case when s — 0. Shining a 
pulse on the QD we obtain due to the Pauli blocking ef- 
fect in QDs the usual two-level situation: no fluorescence 
from initial state |0), full fluorescence from state |I). The 
final state measurement i.e. measuring the spin state of 
the excess electron in the QD is obtained by the quan- 
tum jump technique (e.g. Ref. |I4|'): when the original 
state of the spin in a QD is |1) a fluorescence pattern is 
obtained, whereas the state |0) is completely decoupled 
from the laser field since exciton creation is blocked by 
the Pauli principle. 

The typical time scale which limits the process is the 
spin coherence time in the QD which is of the order of 
microseconds, i.e. in a time of that order of a microsec- 
ond the spin of the electron in the QD will flip from the 
|1) fluorescing state to the |0) dark state and the fluores- 
cence pattern will be terminated. The average number of 
photons emitted before the spin typically flips its state is 
given by the ratio of the spin coherence time to the typ- 
ical rate for spontaneous emission, which is of the order 
of a few nano-seconds. Therefore, typically one should 
obtain on the order of 10'^ photons before the original 
spin information is destroyed. 



B. Case with mixing 

A realistic QD will exhibit mixing of the heavy and 
light hole states. This invalidates the assumption of per- 
fect Pauli blocking with a-~^ light and can be viewed as 
a rotation by an angle —e in the {|0),|1)} space. The 
mixing parameter e will typically be of the order of the 
lattice constant, a, to a typical length scale defining the 
QD in our case a/L « 0.1 where L is the size of the dot 
in the growth direction. 

Introducing mixing requires one to treat the full three- 
level lambda configuration shown in Fig. ^1 As opposed 
to the usual atomic lambda configuration ||4^], here one 
can not distinguish between the |0)(a::| and |l)(a;| tran- 
sitions. These two transitions are mediated through the 
same photon. The dissipative evolution of the density 
matrix p{t) is given by |4j| 



Pxx = iiT [Pix - Pxi+ £ (pox - Pxo)] - (1 + £^)kPx 
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(81) 



K 



The probability that at time t no photon has been emit- 
ted, starting from state a at time to, is 



P^^\t-to)^tr [p{a,t)]. 



(82) 



where at the initial time we take p{a,tQ) = |a)(Q;|. 
Fig. II II shows an example of their evaluation with f2 = 3 
meV, K = 1 ns^^ and e = 0.1. 

Note that, in contrast to the "common" lambda con- 
figuration in Eq. both of the recychng terms, 
Kpxx and e'^Kpxx, are missing, since it is the same pho- 
ton that induces both these transitions, i.e. we can not 
distinguish between the two transitions via photon detec- 
tion. This implies that, when the first photon is emitted, 
say at time ti, the system collapses cither into state |0) 
- with probability po — e^/(l -I- e^) - or into state |1) - 
with probability pi = 1/(1 -I- e^) -, whence the evolution 
starts over again. Therefore the probability that, at the 
time t > ti {i > 1), the (i + l)-th photon has not been 
emitted, is 



U) 



1 



(83) 



which is independent of the initial state \a). A typical 
photoemission pattern will look like Fig. El a sequence 
of pulses, each one made out of a bunch of the order of 
1/e^ photons, separated by no-emission windows. This 
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FIG. 11; Probability that at time t the first photon has not yet 
been emitted, starting from state |0) (above) or |1) (below) 
at time t = 0. Parameters are quoted in the text. 



is the typical quantum- jump pattern one obtains in the 
presence of an emission probability having the form of a 
sum of different exponentials like Eq. H83|) . 
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and heavy hole states as described in Sect. lVIl To be spe- 
cific, let us refer to the left part of Fig.|Sl In that case, the 
transition to be excited for probing the QD state is the 
one marked as (2). Hole mixing results in an unwanted 
coupling to the transition (1). This can be compensated 
for by simply adding a small component of ct_|_ light, pro- 
portional to the mixing parameter e. The error affecting 
this operation would scale linearly with the imprecision 
in e, which is not straightforward to predict theoretically, 
yet it can be measured in a real situation to a good ac- 
curacy. 



C. Case of perfect detection 

We start by considering the error for the case r] = \ 
where t] is the parameter describing the detection effi- 
ciency. On introducing hole mixing the detector would 
still be able to discriminate between the two logical states 
but mixing will be the cause for two types of measure- 
ment errors which can occur. Starting with the system in 
state |1) there is a possibility for no photon to be emitted 
from the QD during the whole measurement time. The 

probability for this type of error is given by: P[°\t). In 
the other case starting with the system initially in |0) at 
least one photon might be emitted during the measure- 
ment time. The probability for this sort of error is given 
by: 1 — P^\t). The measurement time has to be cho- 
sen in such a way as to minimize the sum of these two 
errors. For the same parameters employed in Fig. 1111 we 
obtain an estimate for the optimal measurement time of 
the order of a few tens of ns. What typically happens in 
practice is that, as shown in Fig. 1121 by appropriate time 
windowing the first bunch of photons coming from state 
|1) can be safely discriminated from the (later) photons 
coming from state |0). 



I 

50 100 150 

t (ns) 

FIG. 12: Simulation of photon counts for a system starting 
from state |0) (top) and from state |1) (middle). An expanded 
view of the first few photon counts is displayed in the bottom 
graph. Parameters are the same as in Fig. 1111 

The only feature which allows for discriminating the 
two patterns is the first bunch of photons, which arc 
emitted almost immediately in the case of state |1), and 
after a sensible delay in the case of state |0), due to the 
fact that, prior to the first photoemission, it was still 
P^^\t) ^ P[^\t). Therefore, a detector with 100% effi- 
ciency would be still capable of discriminating between 
the two logical states even in the presence of hole mixing. 

Another option would be available in II- VI semicon- 
ductor systems, showing energetic inversion between light 



D. Finite detection efficiency 

We now consider the case in which 77 < 1. The lowest 
detector efficiency in which we can still hope to discrimi- 
nate between the two logical states is given by 77 = 1/{N), 
where [H (N) = l/e^ is the average number of photons 
to be emitted before a system starting off in state |1) 
fiips to state |0). In our case this is not a tight con- 
straint, since semiconductor photo-detectors have a very 
high quantum efficiency ijmat ~ 0.98. The typical 
wave length emitted by the recombination process in the 
QDs lies well within the spectral window which is due to 
the cutoff by band gap energy of such detectors. The 
main source for low detection efficiency is due to the 
probability for the emitted photon to reach the detec- 
tor, i.e. the difficulty arising due to finite angle coverage 
of the detector. The situation however can be signifi- 
cantly improved by coupling the QD with a microcavity 
as described in (14| . 
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It is important to note that we can use an avalanche 
photon counting mode so that each photon arriving cre- 
ates one e-h pair which then amphfies in the device to 
produce a current spike. The need to wait a few nanosec- 
onds before detecting the next photon is not a hmitation 
in our case since the measurement process we are consid- 
ering is essentially a one shot measurement process, as 
long as the dark count is low enough. 

Working with a detector with a finite efficiency means 
that we have to choose the measurement time so as to 
ensure the fluorescent state emits a few photons thus in- 
creasing the probability one of them will be detected. 
This increases the probability for an error due to a photon 
being emitted by the initial state |0) since Plf\t) decays 
exponentially with time. Moreover a further possibil- 
ity for error is introduced into our measurement scheme. 
Starting initially in state |1) the QD can emit a pho- 
ton/photons which will go undetected and the spin can 
flip into state |0), i.e., the information regarding the spin 
state is lost without being detected. The probability Pe 
for such an error is given by 



realistic conditions and evaluated the corresponding fi- 
delity. Our scheme turns out to be able to suppress the 
effect of both of these decoherence sources on the pro- 
posed gate, and therefore constitutes a viable proposal 
for all-optical quantum information processing in semi- 
conductor quantum dots. 
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APPENDIX A: COUPLING CONSTANTS: 
MICROSCOPIC CALCULATION 
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(84) 



which is simply the sum over n incidents in which the 
emitted photons were not detected and no spin-flip oc- 
curred and on the n-\-l incident such a spin flip occurred 
(without the photon being detected). This type of error 
turns out not to be particularly sensitive to the time of 
measurement. Given that the number of photon emitted 
in the flrst bunch in Fig. El is of the order of e^^ w 10^ 
in our case as discussed above, and taking an efficiency 
r] = 0.9 we obtain an error due to flnite detection effi- 
ciency of the order of 0.1%. 



IX. CONCLUSIONS 

To claim that a certain implementation scheme for 
quantum information processing is viable, one has to 
carefully understand the fundamental sources of deco- 
herence acting in that specific physical system, and to 
show that they can actually be controlled. To this aim, 
in this paper we analyzed in detail the different decoher- 
ence mechanisms affecting a recently proposed all-optical 
scheme for quantum computation based on electron spin 
in quantum dots. In particular, we took into account the 
effect of hole mixing and of coupling to phonons at a fi- 
nite temperature, estimating their impact on each of the 
building blocks of a quantum computer: single- and two- 
qubit gates, and state read-out. We developed a strat- 
egy to circumvent such unwanted effects via an adiabatic 
laser excitation scheme, simulated its performance under 



The phonons are coupled to the charge distribution in 
a QD by means of either deformation or piezoelectric cou- 
pling potentials. The corresponding coupling constants 
are either 



qD{q) 



V2pw(q)F 



in the case of deformation coupling, or 



(Al) 



(A2) 



in the case of piezoelectric coupling. In both cases p is 
the mass density of the sample, and 



Tifi \ 247reei4 .rQxqyqz 
M{q) = ^qxqyqz = M 1— 



(A3) 



is the coupling potential. Here e is the dielectric constant 
of the sample and eu is the material constant. The form- 
factors 



p(q) = j dv[\Mv)\^ - |Vc(r)n exp(-iqr), (A4) 
and H 

D{ci)^ y"dr[Z?,|V.(r)p-Z?,|^,(r)|2]exp(-iqr) (A5) 

are related to the exciton charge density. The wavefunc- 
tions V'^jS'iid ■0c describe the hole and the electron states 
making up the exciton. Dc, are the deformation cou- 
pling potentials. 

The calculation of the coupling constants relies on a 
microscopic model. We consider a QD in a static external 
electric field Fq directed along the x-axis. Within the 
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simplest model with harmonic confinement potential the 
Hamiltonian for the single-particle electron (z = e) and 
the holes {i — h) states is given by 



H = 



P 

2™,- 



2 2 

mMfrf 



Oi 



(A6) 



where r^i — eiFt^/niibjf is a measure of the electric field 
strength in "oscillator units of length" , uji is the fre- 
quency of the confining potential, and are the mass 
and the charge of the particles (the electrons or holes). 

To clarify the effects of external electric field and com- 
pare the relative strength of the different types of the 
coupling, let us consider first the somewhat unrealistic 
but otherwise simple model of a spherically symmetric 
QD. In this case the wavefunction is given by 



1 



3/2 



exp 



(r - rp, 
2/? 



(A7) 



where U = (miUJi)^^/'^ is the ground state localiza- 
tion length. Then, in the case of deformation coupling, 
a simple calculation gives the following expression for 
J— function (for simplicity we put le = lu)'- 



J{lo) 



UJ 



Dl 



2D,.D, 



sin(2u;ro/u) 



.(A8) 



The obtained result shows the two important features. 
First of all, since ^ Dh the J— function is super- 
ohmic in the external field of any strength. Essentially 
this means that a static electric field does not qualita- 
tively change the interaction with the phonons. Sec- 
ondly, the function J(w) has a large frequency cutoff at 
Lui ^ u/l^ which is nothing else but the inverse flight time 
of a phonon through the QD. 

The piezoelectric coupling behaves somewhat differ- 
ent. Recalculating the J— function using the coupling 
constant l|A2|l we find 



/(2-ro 

u 



(A9) 



where /(O) = 1 and f{x) ^ at x ^ oo. Since the piezo- 
electric potential M is the same for the electrons and for 
the holes, in the limit of small electric field strength we 
have: 



67207rVu7 



(AlO) 



The coupling is still superohmic, but it contains a larger 
power of u> than the deformation coupling. Moreover, 
the value of the coupling potential M is also numerically 
small for common materials and thus the interaction of 



a QD with the phonons is dominated by the deformation 
coupling. 

The large electric field limit of Eq. (|A9|) is somewhat 
interesting. One can see that we have 



2807r2pM3 



(All) 



which is, at first glance, an indication of ohmic type of 
coupling. Nevertheless, one should keep in mind that the 
J^*— function decreases quickly when lu }^ lui and hence 
the argument of the /-function never exceeds ~ ro/l. 
This means that the function f 1 practically every- 
where in the course of the integration in Eq. (|48() and 
hence the coupling remains superohmic in the whole rel- 
evant parameters range. 

The analysis above paves the way to a more realistic 
calculation. Consider the Hamiltonian l)A6|l acting in 
2D [x- and y- directions), whereas the motion in the z- 
direction is confined within a box of length L^- The 
wavefunction (of the ground state) is 



^'^^/l:^^"z:V^M — 2r- 



}, (A12) 



where U — ^Jm^Lol. The Fourier component is given by 

y" d^r^a exp{-iqr} = ^^,(g,)F(q), (A13) 
where q = (g^;, qy) is the vector in the xy— plane, and 

-,2;2 



and 



i^(q) = exp|-^ + «qro}, 
47r^i[exp(ig2L) — 1] 



(AM) 



(A15) 



We note that, in spite of its ugly appearance, the form- 
factor Fz is nowhere singular on the real axis and quickly 
decays when q^L 3> 1. The normalization ensures that 
FM = 1. 

For simplicity consider the limit of very strong con- 
finement: Lz ^ lev Then, neglecting the piezoelectric 
coupling and using the zero-argument value for the func- 
tion Fz , we obtain 



where we have 



^"4-2^1 
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(A16)" 
(A17) 



At small value of its argument this function gives 1 and 
vanishes when x oo. The presented J— function is al- 
ways superohmic and is not qualitatively different from 



20 



that considered above for an idealistic spherically sym- 
metric QD. It shares all the important features of the 
simpler model above. In particular, the large-frequency 
cutoff is defined by the same inverse phonon flight-time 
through the QD: uji ~ u/l. 

In short, we presented a number of examples of the 
spectral function calculations. Both in the case of piezo- 
electric and of deformation coupling, the J-function is 
superohmic, with the exponents s = 3 and s = 5, respec- 
tively. The quantity uji ~ u/l plays the role of the high 
frequency cutoff. 



APPENDIX B: ADIABATIC EFFECTIVE 
HAMILTONIAN: APPLICABILITY OF 
PERTURBATION EXPANSIONS 

The phonon interaction terms in the Hamiltonian H47|) 
may well be large and hence the perturbation theory ex- 
pansion in powers of Aq may not always be well justified. 
A remarkable opportunity to extract non-perturbative 
results originates from our adiabatic assumption. 

Indeed our quantum gate proposal relies on adiabatic 
manipulations of the external parameters of the Hamil- 
tonian (|47|l . Assume that both n{t) and A{t) change on 
a time scale r. Extreme adiabatic condition implies 



OJlT > 1, 



(Bl) 



i.e. the considered QD dynamics can be considered slow 
for almost all of the phonon modes. According to [33 |. 
this condition can be formally used to adiabatically elim- 
inate all the phonon modes with frequencies exceeding a 
certain cutoff lu^,: <^ uj^, <^ uji and obtain the effec- 
tive adiabatic Hamiltonian for the slow phonon modes 
with frequencies cj ^ . At zero temperature T = the 
effective Hamiltonian takes the form: 



ph 



-A 



E 

q,(^(q)<a 



(|e)(5| + |5>(e|) + 



e e 



(B2) 



E 

q,uj(q)<a 



;(q)6]j6q. 



The Hamiltonian (|B2|I has the same form as the origi- 
nal Hamiltonian (|47ll . The effect of the high frequency 
mode is limited to the renormalization of the frequency 
detuning 



A = A- 



1 



Jiuj) 



OJ 



(B3) 



and Rabi frequency (tunneling term in the spin boson 
model) 



= J7 exp 



duj J{lu) 



(B4) 



In the superohmic case the integral in the exponent of 
Eq. HB4p converges at ^ and hence does not de- 
pend on Lu^. Thus, in the adiabatic approximation we 



can conveniently set = 0. This means that at zero 
temperature even a strong phonon coupling leads only to 
the renormalization of the Hamiltonian parameters (the 
Rabi frequency). The measure of the Rabi frequency 
renormalization yields the perturbation theory expansion 
parameter: 



n 



duj J {to) J{ijJi) 



< 1. 



(B5) 



The effective Hamiltonian (jB2|l does not rely on this con- 
dition, whereas a perturbation theory for arbitrary fast 
processes would do. 

The situation is somewhat different at finite tempera- 
tures. The QD can still change its state "coherently" , i.e. 
in the course of Rabi oscillations which are characterized 
by the new temperature-dependent renormalized value, 
called Huang- Rhys factor in (34i |: 



£7 = 51 exp|— j 



duj J [to) 



1 



2iV(c.)]}, (B6) 



where N{x) — [exp(x/r) — 1]^^ is the Bose occupation 
number. This is nothing else but the generalization of 
Eq. (|B4|) now taking into account finite occupation of the 
phonon modes. In addition to that, the QD can change 
its state by either absorbing or emitting a phonon. This 
process is called incoherent tunnelling and has no effec- 
tive Hamiltonian form (generally speaking there is no 
obvious separation between the fast and slow variables). 

In the quantum gate proposal we suggest to operate 
our qubit starting from the ground state and adiabat- 
ically changing parameters. This means that the QD 
cannot emit a phonon of a high energy w > (since it 
is in the ground state) . The absorption of a phonon from 
the thermal bath requires finite occupancy of a state with 
energy lo ^ U,, which can be made exponentially small, 
provided T < fi. 

The results of this section also make sense if com- 
pared with perturbation theory. Consider the case of the 
phonon coupling. There is no corrections to the eigenen- 
ergies (|37|l to first order in powers of A^. The second 
order perturbation theory gives 



XI j[l±COs(0)]2 



VA2 + f72 



(B7) 

In the adiabatic limit (|Bip we can expand in powers of 
A/uji and Q/loi to obtain the expression: 



5E± = i(l±cos6i) J dio 



J{lo) 



sm 



Va2 -hfl2 j 



J{co) 



(B8) 



The same expression could be obtained by substituting 
the renormalized values (|B3IB4ll into the adiabatic ener- 
gies (|37|l and expanding the obtained expession in pow- 
ers of the small parameter (jB5|l . This once again shows 
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that the results of the adiabatic renormaUzation coincide 
with perturbation theory whenever both approaches are 
equally valid. 

We conclude that the effective Hamiltonian (|B2|I with 
renormalized value of the Rabi transition amplitude ljB6|) 
can be used as a good non-perturbative tool to study the 
QD dynamics in the presence of phonons. Its validity is 



ensured by the fact that the J-function for the deforma- 
tion coupling is superohmic and relies on the timescale 
separation condtion IjBip . In the case when ljit > 1, the 
adiabatic approximation fails and one has to resort to 
perturbation theory, whose expansion parameter is given 
by 
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